NASA Technical Memorandum 102078 
AIAA-89-2906 


Transient Flow Thrust Prediction for 
an Ejector Propulsion Concept 


Colin K. Drummond 
Lewis Research Center 
Cleveland, Ohio 


(NASA-TM-102078) TRANSIENT FL0H THRUST 
PREDICTION FOR AN EJECTOR PROPULSION CONCEPT 
(NASA. y Levis Research Center) 10 pCSCL 21E 


G3/07 


Prepared for the 

25th Joint Propulsion Conference 

cosponsored by the AIAA, ASME, SAE, and ASEE 

Monterey, California, July 10-12, 1989 


N89-243 18 

Unci as 
0210776 




Transient Flow Thrust Prediction for an Ejector Propulsion Concept 


Colin K. Drummond* 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


Abstract 

A method for predicting transient thrust augmenting 
ejector characteristics is introduced. The analysis blends 
classic self-similar turbulent jet descriptions with a mixing 
region control volume analysis to predict transient effects 
in a new way. Details of the theoretical foundation, the 
solution algorithm, and sample calculations are given. 

Nomenclature 

A cross-sectional area 

b jet half-width 

B channel width 

C i coefficient defined by equation(24) 

KE kinetic energy 

m mass flowrate 

m cy control volume mass 

M momentum flux 

M cv control volume momentum 

M ] p Mach number of primary nozzle flow 
P pressure 

t time 

v velocity 

x transverse coordinate 

2: streamwise coordinate 

Z, integral of defined function 

p density 

function defined by equation(4) 
a turbulent flow constant, equation(20) 

£ dimensionless transverse coordinate 

subscripts 

cv control volume 

e entrained stream 

m primary stream centerline 

1 s secondary steam at station 1 

1 p primary stream at station 1 

Background 

Research on design methodologies for integrated air- 
craft and propulsion flight control systems requires accu- 
rate subsystem component simulations. In principle, these 
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simulations must mimic steady-state and transient 
component effects. A NASA Lewis research program is 
currently underway to develop a "real-time" simulation, 
including system transients, for short take-off vertical 
landing (STOVL) aircraft. The thrust augmenting ejector 
shown in Fig. 1 is considered a potentially valuable pro- 
pulsion subsystem element for the powered-lift aspect of 
STOVL aircraft. To explore ejector concepts further, the 
initial STOVL system simulation requires an ejector sub- 
system simulation. 

For analysis, a thrust augmenting ejector is usually 
divided into a nozzle, inlet, diffuser and mixing regions. The 
ejector shown in Fig. 1 is a mechanically simple fluidic pump 
composed essentially of two components: (1) a "primary" 
jet nozzle issuing into (2) a shroud. This arrangement 
permits entrainment and acceleration of a secondary flow 
(within the shroud) by the primary jet. A diffuser section 
attached to the shroud allows control over the ejector 
discharge conditions. 

An ejector simulation that includes ejector transients 
and with the potential to run real-time is not currently 
available; the motivation for the present work is the need 
to develop one. Of particular interest is the mathematical 
treatment of the ejector mixing region. 

Motivation for a New Modelling A pproach 

Motivation for a new perspective on thrust augmenting 
ejector flows extends from the failure of existing mathe- 
matical models to meet current integrated flight / pro- 
pulsion controls (IFPC) simulation requirements. A 
successful approach must satisfy three key constraints: 

1. The formulation must handle unsteady flows and 
include a model for turbulent flow 

2. The mathematical model must be predictive (not 
parametric) in nature, 

3. The final system of equations must be amicable to a 
real-time simulation objective. 

The first requirement points out the need for some type of 
characterization of turbulent mixing and entrainment 
phenomena in the mixing region; a description of these 
internal effects is essential for transient-type analyses. To 
meet the second requirement, only two data sets should be 
prescribed: (a) the primary jet "control-valve" setting (with 
associated thermodynamic data), and (b) free-stream 
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atmospheric properties. From this, the secondary flow inlet 
conditions and the ejector thrust augmentation (as a 
function of time) are predicted. The third requirement 
emphasizes the ejector simulation must not become a 
bottleneck in the STOVL simulation (and thereby elimi- 
nating many traditional CFD methodologies in this appli- 
cation). 

Constructing a mathematical model within these con- 
straints is not a simple task; existing simulation approaches 
accommodating this scope were not found in the open 
literature. Traditional methods of analysis are described by 
Porter and SquiresC 1 ) and Addy and Dutton(^). 

Problem Solution Overview 

It is clear the IFPC ejector simulation task demands 
innovative approaches to unsteady turbulent flow treat- 
ment. Research at NASA Lewis has yielded a method for 
ejector analysis that represents a potential solution to the 
simulation problem. An empirically based model for the jet 
mixing turbulent interaction region is explored within the 
framework of a control volume analysis. This approach 
provides a rational foundation for the introduction of 
steady-flow data to "calibrate" an unsteady flow simulation. 
Preliminary performance predictions are intuitively rea- 
sonable, though verification of accuracy awaits availability 
of transient ejector performance data. 

Here, a transient flow is a ’temporary’ unsteady flow, 
associated with, for example, a change in ejector operation 
from one steady-state condition to another. Contrast this 
with oscillatory flows in which a periodic time-asymptotic 
flow character is exhibited. Ejectors utilizing pulsed pri- 
mary nozzle flows are of the latter type. In the present 
simulation the focus on transient, not oscillatory, ejector 
phenomena descends from flight-critical aircraft flight 
control scenarios; an example would be transition to for- 
ward flight from vertical take-off. 

Mixing Region Focus For the analysis of a given 
unsteady flow problem, it is quite convenient mathemat- 
ically if a quasi-steady formulation can be assumed valid. 
This leads to the use of steady flow equations in an unsteady 
flow analysis; at each instant in time the flow is assumed to 
instantly respond to changes in the boundary conditions. 
The validity of such an assumption, however, depends on 
whether the characteristic time of the forcing function is 
the same order of magnitude as the relaxation time of the 
flow. In the present work only the inlet and diffuser benefit 
from the quasi-steady flow assumption; this directs our 
focus more appropriately on the control volume analysis of 
the mixing region. The notion here is that the physics of the 
mixing region are the cause of a situation that yields the 
pressure gradient effects in the inlet and diffuser. 

Mode! Overview The mathematical framework of the 
present model is essentially a blend of three basic ideas. 
First, identify a finite number of control volumes (about 
five) of fixed size, partitioned only in the streamwise 
direction (the size of the control volumes are time inde- 
pendent). Translation of the Navier-Stokes equation into 
finite volume form provides three times as many algebraic 


equations as finite volumes. The robustness of the finite 
volume approximation minimizes the number of finite 
volumes needed; this is synonymous with reduced execution 
time. 

Second, introduce the assumption that dimensionless 
self-similar velocity, temperature, and pressure distrib- 
utions are applicable to the primary ejector flow and, 
furthermore, that the classic Abramovich^) self-similar 
turbulent flow expressions are valid. This step provides the 
simplest turbulent model applicable to ejector flows. As a 
first approximation to meeting the primary and secondary 
jet mixing interface condition, the static pressure is assumed 
uniform in the transverse direction, but not in the stream- 
wise direction. Characteristic profiles (see Fig. 2) provide 
a convenient link between the primary and secondary flows 
within the mixing region elements; Fig. 3 illustrates that 
there now results a very simple "virtual grid" on which to 
base the finite volume computations. 

The third idea to introduce is that the kinetic energy 
gain of the secondary flow is expressible as a function of the 
kinetic energy loss of the primary flow. Contrast this with 
the traditional approach (see Korst and Chow^ ) of 
expressing the total kinetic energy of the secondary flow as 
a function of the total kinetic energy of the primary flow. A 
new kinetic energy function is needed since the traditional 
approach provides, by default, a quasi-steady flow result. 
In both cases an empirical constant must be introduced; in 
the traditional approach the constant derives from obser- 
vation of normalized steady flow data — in the present work 
the constant is determined by matching asymptotic thrust 
predictions for a given ejector performance window. 

In summary, balances of mass, momentum, and energy 
are laminated on a simple virtual grid to provide the 
governing equations of the problem. Unlike many control 
volume analyses, the unsteady terms (time rate-of-change 
of, for instance, momentum of all particles inside the control 
volume) are not dropped in the present approach so we are 
forced to estimate field conditions in the ejector mixing 
region. In this work, the self-similar turbulent profiles 
approximate the velocity distribution within the control 
volume and the balance of energy divided into the (usual) 
thermal energy balance and a new function for mechanical 
energy. 

Formulation 

Mixing Re gion Control- Volume Equations 

A traditional one-dimensional flow approximation 
yields the following mass and momentum conservation 
laws: 

Mass 

= Y'Tht (1) 

Momentum 

“fp vdV - - pv(v • n)dA - <£ n • (pI)dA (2) 
dtJv ~ J A - - — J A~ 
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From a philosophical point of view, the presence of 
surface integrals (that in steady flow analyses yield sums of 
average quantities) reflects an increase in surface flux 
complexity, but allows for simplification of the volume 
integral representation. Control volume equations for 
energy and entropy balances are given by Thompson^); for 
brevity of presentation the description of the finite volume 
method of this work will be restricted to eauations(l) and 
(2). More details are given by Drummond W, 

Turbulsatist Approximation 

The present turbulent jet geometric approximation is 
| shown in Fig. 2. Extending downstream from the mixing 

region inlet plane there exists a potential-core region 
i characterized by a fairly uniform centerline velocity, with 

no transverse component. This is distinguished from the 
mixed-flow region where the centerline velocity decay arises 
from momentum transport to the entrained fluid. 

Self-Similar Profiles In the characterization of an 
element of the mixing region, the velocity assumes a 2-D 
planar turbulent jet self-similar profile of Abramovich^) 
for co-flowing jets: 

u = u c ( 1 - 4>) + - /(*) (3) 


m, - 2I/MX0.45u m + 0.S5u, + (S-lK)< 

- 2Vb i Z l (8) 

The time-derivative term in the mass conservation equation 
is 

(£i - — (!?), <*> 

where the characteristic density for the finite volume is now 
approximated by the value of the density at station j (in 
practice this is a good assumption as long as field variable 
gradients are 'modest" in size). If the characteristic jet 
expansion width, b, also assumes its value at j, then 
substitution and re-arrangement of the continuity equation 
yields 

(^) - * ~ tilbi (10) 

\dt), BAz ' J 

Computation of the jet half-width, b, derives from the 
momentum equation for incompressible flow applied to 
the first finite volume; a rectilinear jet expansion is assumed 
therefrom and compressible flow restored for subsequent 
calculations. 


where 

*(?) 


/(l-? I5 ) 2 ,0<^<l\ 

\ 0, !<?<i / 


Application to the Momentum Equation 

When the self-similar profiles are substituted into 
( 4 ) equation (2) and the algebraic dust settles, the jet centerline 

velocity derivative is 


and 





M k v m + v„F l 

(dp) 

= x/b 

(5a) 

l dt Jj 

2VbjAz p 

[ dt J y 


| = B/b (56) 

Inclusion of the core region requires a slight modification 
of these profiles; details are given in Abramovich^). 

Since there is a static pressure matching condition at 
the jet boundary, a uniform transverse pressure distribution 
is assumed; in the longitudinal (axial) direction a finite 
pressure gradient exists. This feature partially distinguishes 
free-jet and confined-jet analyses. As a first approximation, 
introduction of the uniform pressure profile into the ideal 
gas law suggests the transverse density gradient can also be 
uniform. 


Applic a tion t o C onservation o f Mass 


For the generic sub-region k , bounded by surfaces at i 
and /, mass conservation yields 



( 6 ) 


First, the mass flux is written 

rn< * f PiVidA *21 vf (7) 

J A J 0 

From which integration over the self-similar profiles yields 


where 

M k 


1 + 0.45 
0.45 

2I/6 ( p ( (Z 2 . ( + |F ( ) 


( 12 ) 


-2l/b yPy (Z 2/ + |F y ) (13) 

“ J (14) 

J 0 

Note in equation (11) that the centerline velocity derivative 
is a function of the density derivative, equation (10). 

Although application of the finite volume approxima- 
tion has not been shown here for energy and entropy, the 
procedure is the same as described above. The result is a 
system of three equations describing the time derivatives 
of the jet centerline density, velocity, and temperature in 
terms of the state of the entrained flow. Integration of the 
system commences when the entrained flow state is defined. 


Entrained Velocity Augmentation 

Analysis of the primary and secondary flow interaction 
has not, to this point, been completed. By themselves, the 
self-similar profiles close the loop for steady-state flows , but 
not transient ones ! This section provides an approximation 
for the turbulent flow kinetic energy exchange mechanism 
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to characterize the influence of primary flow changes on 
the secondary flow. When the change in state of the 
secondary flow is established the corresponding secondary 
field variables can then be computed. 

Secondary Stream Energy Chang e Computation of the 
gain in secondary flow kinetic energy due to mixing is given 
by 

A KE = (15) 


of its total loss. In balance, however, the total change of 
kinetic energy of the primary flow is greater than that of 
the secondary flow. 

In the works of Korst and ChowW and Chow and 
AddyO) the relationship between the change in entrained 
flow kinetic energy and the total primary flow kinetic energy 
is given by 


where ^defines the jet boundary streamline (for which the 
primary mass flow through station / is equal to primary mass 
flow through ;). For the present discussion this dividing 
streamline position is assumed known; Korst and ChowW 
discuss the typical approach of analysis. Expanding the 
equation for the change in kinetic energy yields 

A KE = (v 3 -vv 2 ,)d% + j\v 3 -vv 2 )d^j 

(16) 

Substitution of the self-similar profiles into this expression 
and integrating the result provides 

A KE = N\Sbp{v 3 9 (H ^ + H^- F 3 - F 4 ) 

+ v 2 m v 9 H 2 +vlH 4 } (17) 

where the detailed expressions for the integrals H\ and F\ 
are given by Drummond(^). It should be noted that in the 
present work the integrals are independent of time. 

Primary Flow Energy Change Similar to the way in 
which the change in secondary flow kinetic energy was 
computed, the energy loss of the primary flow is given by 

AKE - J\ US) 

where the limits of integration reflect interest in the domain 
of the primary jet cross-section. 

Evaluation of the integral at station i yields the result 

A KE = Vbp {v\H x + v 9 v\H 2 

+ v* m H< 

-v 3 m F,-vy m F ,) (19) 

Again, reference(S) provides the expressions for the inte- 
grals implied in Fj and H[ 

Kinetic Energy Balance Computations for a specified 
steady-state condition show that the change in kinetic 
energy due to mixing is not the same for the secondary flow 
as it is for the primary. In fact, the gain in kinetic energy of 
the secondary flow is entirely due to the mixing process, 
while the mixing loss of the primary flow is only a fraction 


where a value of 12 for a for turbulent flow traditionally 
provides a reasonable match between theory and experi- 
ment. A more accurate relationship includes the effect of 
the primary flow Mach number, 

o - 12(1+ 0.23M ) (21) 

The important feature of this result is that the change in 
secondary flow kinetic energy has the functional form 

A KE ls = (22) 

The difficulty with the energy transfer function 

described above is that it provides (by default) a quasi- 
steady flow approximation. It therefore cannot be used for 
the transient flow in its present form. To entertain local 
changes in primary-secondary kinetic energy in a way that 
does not burden the computational procedure, consider the 
modified function 

A KE ls - (23) 

where the subscript m denotes the change in kinetic energy 
due to mixing. The proposed form of this function is given 
by 

A KE\T’ *KE\'“ m * C , ( A KE\,' m (24) 

This function for secondary flow KE enhancement arises 
from the assumption that local velocity gradients in the 
steady flow case are typically less than the gradients 
experienced in the transient mode. Here, the introduction 
of an engineering approximation also results in the intro- 
duction of an undetermined constant, Cj . The alternative 
is to establish N computations of the kinetic energy 
exchange to coincide with the N control volumes of the 
mixing region; the present method permits post-processing 
KE information at the completion of mixing region calcu- 
lations. Numerical experimentation indicates 

0. 1 < C j < 0.5 (25) 

In summary, knowledge of the total change in secondary 
flow kinetic energy permits an update to the secondary flow 
state that corresponds with the dynamics of the primary 
flow. Once the velocity of the secondary flow has been 
extracted, the new inlet density and pressure are computed 
from an ideal flow inlet analysis. 
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Numerical Method 

Transformation of the control volume equations into 
finite-volume form provides a set of equations for the 
time-derivative of the field variables for each mixing region 
control- volume. These results participate in the solution as 
follows: 

1. Define the computational grid (5 or so elements), 

2. Initialize field variables by setting time derivatives to 
zero for mass, momentum, and energy, then solve for 
the associated (steady-state) field variable distribution, 

3. Compute field variable time derivatives for each mixing 
region finite volume (use field variable values, as 
needed, from the previous time step), 

4. Advance field variables forward in time explicitly using 
first order finite differences for the time-derivatives, 

5. Update the entrained velocity, 

6. Repeat steps 3-5 until desired time is reached. 

Figure 4 summarizes this procedure. 

Results 

Characteris tic Test Case 

In the absence of data from transient flow ejector 
experiments (or even from modern multi-dimensional 
Navier-Stokes solvers), "verification" that the proposed 
ejector analysis can provide reasonable thrust predictions 
must be left to engineering judgement. Because of this, a 
"familiar" ejector forcing function must be used. In the 
present work the system response to a step-function input 
is not only a characteristic transient case study, but the 
scenario is also coincident with typical STOVL ejector 
application. For demonstration purposes the ejector system 
response to a step-change in primary nozzle flowrate from 
18.7 to 21.85 lb m /sec is chosen because (a) experimental 
steady-state data at each of these operating points is 
available, and (b) the 17% change in primary flowrate is 
well beyond a "small"-perturbation examination (this 
exercises the system non-linearities). 

Calculation Results 

For the mixing region finite-volume length of 0.18 ft and 
a characteristic mixing region velocity of 500 ft/s, the 
characteristic time step for computations is approximately 

0.4 ms. To avoid infringing on this stability limit a compu- 
tational time step of 0.1 ms was chosen; 100 time steps 
provided the necessary interval for examination of the 
step-function test case. 

The empirical coefficient in the transient analysis (re- 
quired for calibration of the primary-to-secondary kinetic 
energy exchange mechanism), was selected to match the 
asymptotic transient thrust prediction with the quasi-steady 
value at the new set point. Figure 5 illustrates the predicted 
ejector thrust profile for a coefficient Cj of approximately 
0.25. It appears the 2 millisecond residence time of the flow 
(elapsed time for primary nozzle flow to reach the diffuser 
exit plane) is slightly less than the 3 millisecond interval 


for the thrust to reach a new maximum. Oscillations in 
thrust after that point appear to settle in about 5 millisec- 
onds. 

An unexpected feature of the thrust profile is the dip in 
thrust immediately following the step-change in primary 
nozzle efflux. Examination of the field variable profiles 
reveals this is not a numerical problem, but that the increase 
in static pressure associated with the instantaneous change 
in driving flow temporarily impedes the secondary flow. 
After a short period, the secondary flow kinetic energy 
builds (commensurate with the increase in primary flow 
energy) to overcome this effect, then continues in the 
intuitively expected manner. 

Remarks 

A distinctive second-order flavor is displayed by the 
predicted thrust profile; under a second-order assumption 
the ejector test case has approximately a 0.75 damping ratio 
and a natural frequency on the order of 300 Hz. Although 
the results seem reasonable, it is necessary to conduct more 
extensive computational tests before conclusions about the 
order or linearity (about the perturbation) of the system 
can be made. 

Summary 

An ejector simulation method that includes ejector 
transients and with the potential to run real-time has been 
presented. The finite volume method permits rapid evalu- 
ation of the time dependence of field variables in a thrust 
augmenting ejector mixing region. Accuracy normally 
increases with an increase in the number of finite-volumes, 
though if the intent is to compute thrust, then quite 
reasonable results can be obtained with as few as five finite 
volumesb. 
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Fig. 1 Typical finite volume for ejector mixing region 
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Fig. 2 Nomenclature for finite volume analysis 
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Fig. 3 Computational grid 



Fig. 4 Outline of solution Procedure 
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